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Recent results from PAMELA and ATIC hint that ©(TeV) dark matter (DM) is annihilating, in 
our galactic neighborhood, predominantly to leptons. The annihilation rate is much larger now than 
during freeze-out, one possible explanation of this is a low-velocity enhancement of the annihilation 
cross section. In a model independent fashion, we show that in this case the rate of neutrino emission 
from the Earth, due to DM annihilation, may be greatly enhanced while the rate from the Sun is 
unaltered. There is potential for IceCube to see these earthborn neutrinos while the same parameter 
space will be soon covered by direct detection experiments. Combining these near-future data will 
allow extraction of valuable information about the DM sector dynamics. 



Approximately 20% of the matter-energy budget of the 
universe is due to Dark Matter (DM). The favored can- 
didate for the DM particle is a thermal relic with anni- 
hilation cross section {av) « 3 x 10~^^cm^s~^, a weakly 
interacting DM (WIMP). Many experiments are under- 
way to probe the DM, either directly through its interac- 
tions with Standard Model (SM) particles or indirectly 
through its annihilations to SM particles. Recently sev- 
eral indirect detection experiments have reported results 
which may be interpreted as due to DM annihilations, 
although they could also have an astrophysical origin [l[ . 

The ATIC experiment has reported an jxcess of 
electron-positron flux around 300 — 800 GeV |2j . In ad- 
dition PAMELA 3] is seeing an increase of the positron 
fraction around energies of 10 — 80 GeV and no corre- 
sponding excess in the antiproton fraction [J]. Taken 
together these suggest that 0(1 TcV) DM, annihilating 
preferentially to leptons, is being observed How- 
ever, the annihilation cross section required to explain 
the excesses is substantially larger than 3 x lO^^^cm^s^^ 
[1, Q. The enhancement may be due to a boost factor, a 
nearby clump of DM or a low- velocity Sommerfeld effect 
(see also [l3. El for an alternative explanation) . 

These results are exciting and surprising, not only are 
we possibly observing WIMP DM but maybe also a non- 
trivial DM sector, whose dynamics seems to imply an 
epoch-dependent annihilation cross section. Looking at 
the photon spectrum [H [H, [llJllJli 

and additional 

cosmic ray (CR) experiments 17, 18, 1^ will test this 
emerging paradigm 7] and whether these excesses are ac- 
tually due to DM. However, the photon and CR flux de- 
pends on the DM halo profile and the propagation model. 
We demonstrate here that the same DM sector dynamics 
may induce dramatic changes in the neutrino flux from 
the Earth which give a very different probe of the same 
microscopic phenomena. For a possible signal from galac- 
tic neutrinos see [10, lU . 

Dense bodies such as the Sun and the Earth gravi- 
tationally capture DM particles in their core, resulting 
in a DM density significantly higher than in the galac- 
tic halo. They eventually start to annihilate into SM 
particles, among which neutrinos can escape and travel 



to Earth-based detectors. The resulting flux depends on 
the capture and annihilation cross sections, unless the 
DM has already reached equilibrium which leads to a 
maximal flux, exclusively controlled by the capture rate. 
The effects which enhance the annihilation rate would, 
generically, not affect the capture rate. For instance, an 
ultra light particle with sizable coupling to the nuclei is 
required to Sommerfeld enhance the capture rate which 
is probably in conflict with various precision data. While 
the capture rate remains unaltered, a larger annihilation 
cross section will shorten the typical time for the DM to 
reach equilibrium. Our key observation is that since it is 
very probable that the Earth has not yet reached equilib- 
rium for a relic annihilation cross section [13, 113 , this ef- 
fect would yield order of magnitudes enhancement in the 
neutrino flux from the core of the Earth at IceCube [1^ . 
Moreover this flux will also be correlated with the DM di- 
rect search experiments [251, ^26*1 . A combination of these 
two data sets yields fairly clean information about the 
microscopic nature of the dark sector dynamics. 

Neutrinos via DM annihilation. The competition 
between capture and annihilation of the DM leads to a 
present day DM annihilation rate |27j] 
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where ~ 4.5 x 10^ yrs is the 

A = {(Tv)/Veff, V,ff = 5.7 X 10^2 



age of the Earth, 
e//, Veff — o.i K ru~cm^ (TeV/m^)'^^^ is 
the effective volume of the core of the Earth [1^ and C 
is the capture rate. For the Earth the capture rate is 
dominated by the spin independent (SI) elastic scatter- 
ing (see [1^ for the inelastic case) of the DM off various 
elements in the Earth [27l. 
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where the the sum is over the elements O, Si, Mg, S, 
Fe and Ni, only 3% of the mass of the Earth is ne- 
glected. The DM mass is denoted to^, pg g and tijyo 
are the DM energy density and velocity in the halo in 
units of 0.3 GeV/cm^ and 270 km/s respectively, while 
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the factor fi accounts for the mass fraction and distri- 
bution profile of the element i [l^], whose cross section 



with DM is denoted tr^'. Direct detection experiments 
probe SI cross section of DM off proton, Cgj. To better 
than 1%, (TgJ w -^^'^si ^'-'^ nucleus of mass number 
N. Hence, 



Cm ~ 9.6 X 10"s-i 



(1-270) 
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The maximum rate of DM annihilation occurs after 
equilibrium is reached and is entirely determined by 
the capture rate, Feq — C/2. For times shorter than 
the equilibrium time t^q ~ 1/VCA the abundance is 
grows linearly with time and the annihilation rate is 
Tneq ~ ^AC^t^. With a typical thermal relic annihila- 
tion cross section, ~ 5.3 x 10~"'^s~^ (m^/TeV)^^'^, 
the Earth is far from equilibrium {t^ ^ icq) and not 
a good source of DM-neutrinos. However if the observed 
electrons/positrons excesses are due to a low- velocity en- 
hancement, R, the annihilation cross section can be far 
larger than that of the early universe, A^ = RA^, bring- 
ing the Earth towards equilibrium today. The maxi- 



mal enhancement in the rate is roq/r„ 



which can be several orders of magnitude and is obtained 
for R > {Aj-C^t'^)^^ . The escape velocity at the cen- 
ter of the Earth is approximately 15km s^^ whilst DM 
in the halo has a Maxwell-Boltzmann distribution with 
Wo = 270 km s~^. The Sommerfeld enhancement grows 
as ~ 1/v although this growth saturates at very low ve- 
locities a further increase beyond t; = 270 km s~^ may 
yield more non-trivial information about the DM sector. 
Thus, the enhancement may in fact be even larger than 
that for DM in the halo. It will be useful to define the 
critical capture rate for the Earth: 



= 1/A,tl ~ 9.93 X lO^^g-i f—\ ^ 



(4) 



above which the Earth would already have reached equi- 
librium and boosting the annihilation cross section will 
not result in an enhanced neutrino flux. Direct searches 
experiment such as CDMSII put an upper bound [2^ on 
the SI elastic scattering cross section of 3.5 x 10~ ''pb for 
rriy. = ITeV. Thus, C® < lO-^C^, and the Earth is 
probably still far from equilibrium. 

The capture rate ([3]) is derived assuming that the DM 
velocity distribution as encountered by the earth is Gaus- 
sian. It is possible that in the solar system it differs 
from Gaussian [23, 29] , particularly at the low velocities 
necessary for capture in the Earth and Sun. The DM 
abundance may also differ significantly from the galac- 
tic halo density (see e.g. [23, _2£, ^1, 32] and Refs. 
therein). Both direct and indirect detection experiments 
probe the same nuclear scattering cross section but they 
are sensitive to the different parts of the velocity dis- 
tribution, high and low velocity respectively. Assuming 
a Gaussian distribution allows observations from direct 



and indirect experiments to be straightforwardly corre- 
lated. Furthermore, a future signal at direct detection 
experiments would directly probe velocity distributions 
(through differential energy information) of the DM par- 
ticles [21, m [H at ranges of roughly 40-150kms"^ 
Of particular importance are the Xe based experiments 
which have the lowest threshold, down to approximately 
three times the earth escape velocity 34, 35]. 

Annihilation into primary neutrinos. The muon 
flux at the surface of the Earth is given by: 
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with pp^n the number density of protons and neutrons in 
the medium, respectively 5/9A^ACm^'^ and 4/9iVyicm~^ 
for ice, where iV^ ~ 6 x 10^^ is Avogadro's number. 
da^'" /dEp_ are the weak scattering cross sections of (anti- 
) neutrinos on nucleons 
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where a^-^ = 0.25, 0.15, b^'P = 0.06, 0.04 and a"'^ = bP^''', 
b"'P = a^'" [sH]. The distance R^{E^), the muon range, 
defines the distance traveled by a muon until its en- 
ergy drops below the energy threshold Eth of the de- 
tector, due to losses in the medium. Approximately, 



Rn {En 



with p the density of the 
2.0 MeVcm^g-i 



medium (~ Igcm"'^ for ice) and a 
and (3 ~ 4.2x 10~^cm^g~^ for ice. At IceCube, the energy 
threshold is about 50 GeV and for E^ ~ TeV, the typi- 
cal muon range is a few kilometers, which is longer than 
the detector typical size |48i]. Since the DM is almost at 
rest, the muon neutrino flux at the surface of the Earth 
is monochromatic, d^^^/dE^ = 5{Ei, — m^)Bjy^T / AttR^ , 
with Bf^i, the branching ratio of DM annihilating to neu- 
trino pair and Rq ~ 6.4 x lO'^km, the Earth radius. The 
resulting muon flux is: 
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Combining this with the effective area [37| of the detec- 
tor A^ff{Ef^) gives the event rate in the detector, i.e. 
dN/dEf, = AeffiEf,)d<^>f,/dEf,. This is shown for DM 
masses of 500 GeV and 1 TeV in Figure [1] along with 
the background rate (discussed below) from atmospheric 
neutrinos. 

Secondary neutrino sources. Instead of direct pro- 
duction may be produced from secondary decays of the 
DM annihilation products, and we concentrate here on 
annihilations to charged lepton final states. Muons are 
stopped long before they decay [11] and are not a source 
of high energy neutrinos, whereas taus lose very little 
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FIG. 1: Muon rates from primary neutrinos (red) and at- 
mospheric backgrounds (black) at IceCube. All plots show 
results for a 1 TeV (solid) and 500 GeV (dashed) DM, corre- 
sponding angular cuts have been placed on the background. 
The signal assumes the Earth has reached equilibrium. 



energy and will produce prompt neutrinos. In the case 
where the DM annihilates preferentially to tau leptons, 
which in turn decay into neutrinos, the induced muon 
flux at the Earth surface, taking into account interactions 
with the material in the Earth, can be parametrized by 
the following analytic formula [39|: 



dE„, 
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-mx(p6+P2 cxp(— p3mx)) 



(8) 



where is in GeV. Bfr is the branching ratio of DM 
annihilating into taus and 

« (2x 10-22/W, 0.2, 5x 10-3, 0.1, 6x 10-3, 0.2, lO^^). 

Backgrounds. The main source of background muon 
neutrinos comes from the shower of cosmic-ray inter- 
actions with the atmosphere. The anisotropic induced 
muon flux is then obtained from [30 . ,40j : 
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where 9^ is the zenith angle and the differential fluxes 
of muon neutrinos and antineutrinos are estimated from 
the tables found in Honda et al. [4l|. The function 
R{cose^) = 0.70 - 0.48 cos 6*2 for 9^ > 85° and 1 else- 
where, is the efficiency of IceCube for tracking up-going 
muons. The background can be substantially reduced by 
noting that the signal is coUimated in a cone of half-angle, 
A9 = 1.8°{TeV/E^y/^, about 9^ = 180°, where E^ is the 
energy of the incoming neutrino [13, 112] ■ For the Sun the 
background reduction is limited by the angular resolution 
of IceCube, A9 = 0.5°, which we take about 6^ — 66°. 
While the primary neutrino signal is monochromatic. 
El, = m^, the spectrum of secondary neutrinos is con- 
centrated at low energy due to the hard slowdown of the 




FIG. 2: Statistical significance of a primary (top) and sec- 
ondary (bottom) neutrino signals above the atmospheric 
background. BR denotes Bj^^ (primary) or B^t (secondary). 
Upper curves (red) are the equilibrium fluxes while lower ones 
(black) are the naive fluxes. The blue dot corresponds to the 
critical capture rate, C^, see Eq. @. The vertical lines 
show the present and future upper bounds on (jgj from direct 
detection. The horizontal lines show the discovery reach of 
IceCube. 



DM annihilation products before they decay into neutri- 
nos. Their typical energy is Ei, ^ Eth — 50 GeV, and 
A9 ~ 8° which increases the relevant background by an 
order of magnitude compared to the primary neutrino 
case. Hence, the monochromatic neutrinos offer the best 
hope for a discovery at IceCube. 

Earth and Sun detection potential. The reach 
of IceCube is shown in Figure [2] for both primary and 
secondary neutrinos. We apply energy cuts for primary 
(secondary) signals of 250 GeV < E^ < {Eth < E^ < 
500 GeV). From these plots it is clear that DM that 
does not annihilate directly to neutrinos has very little 
hope of discovery at IceCube, even with a large Sommer- 
feld enhancement, and we concentrate on the primary 
neutrino case. The maximum neutrino flux is given by 
the red line and is well into the 5(t discovery region for 
most of the range that can be probed by direct detection. 
However, since <C icq the expected rate is denoted by 
the black line. Enhancements of > 100 are necessary 
for the ATIC/PAMELA results and niay, depending on 
the details of the resonance structure p, |43| , be consid- 
erably larger for DM in the Earth. Over most of the 
region accessible to direct detection enhancement factors 
of 100-1000 will move us back into the IceCube discovery 
region. Thus, by 2013 we will have probed most of the re- 
gion where neutrinos from the Earth could be discovered. 
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FIG. 3: Statistical significance for primary neutrino from tlie 
Sun as a function of SD scattering off proton. Since the cap- 
ture rate in tlie Sun is more efficient tlian in tlie Eartli it is 
most likely in equilibrium. 



If direct detection experiments make an observation then 
we may have a correlated discovery in IceCube. 

In the core of the Sun the capture rate is dominated 
by the spin dependent (SD) elastic scattering of DM off 
hydrogen nuclei [27| : 

and Q ~ 3.23 x lO^^s-i {TeY/m^f''^ . We em- 
phasize again that there are significant astrophysics un- 
certainties on the DM density and its velocity distribu- 
tions and thus the actual capture rate (we use the 
value from [g^] for concreteness; see also [i^l). Further- 
more, the SD scattering cross section of DM on proton is 
less constrained experimentally, the present bound being 



CT^D ^ 0-8pb (for - ITeV) from KIMS Thus, 
in this case, a future signal from the sun, only from Ice- 
Cube, would be harder to cleanly interpret. It is ex- 
pected, generically, that the SD cross section is 3-4 order 
of magnitude larger than the SI one. Thus, taking the 
above capture rates at face value we see that Cq <C Cq 
for wide range of reasonable DM models. Consequently, 
it is likely that the sun is now in equilibrium and its 
neutrino flux is already maximal, leaving no room for an 
enhancement of the annihilation rate as shown on Fig- 
ure [3l 

Conclusions. Combining the information on the neu- 
trino flux and the direct detection cross section yields a 
fairly robust measurement of the annihilation boost fac- 
tor. The significance of the signal is greatly improved 
in cases where the annihilation channels involve primary 
neutrinos, not an inconvincible possibility [1^. A detailed 
study of a possible signal from primary and secondary 
Sun-born neutrinos may help to determine the primary 
branching ratio. In this case, lack of a related earthborn 
signal would indicate that a low velocity enhancement 
of the annihilation cross section is not the explanation 
for the ATIC/PAMELA excess. Instead, one would look 
for an astrophysics explanation. In the ideal case where 
enough events are observed at IceCube a differential en- 
ergy information could be extracted which may yield fur- 
ther insight towards the DM sector, such as its mass and 
decaying branching ratios. 
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